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Abstract 

We explore a new class of braneworld models in which the scalar curvature 



£T) . of the (induced) brane metric contributes to the brane action. The scalar 

curvature term arises generically on account of one-loop effects induced by 
. matter fields residing on the brane. Spatially flat braneworld models can en- 

o . 

. ter into a regime of accelerated expansion at late times. This is true even 

if the brane tension and the bulk cosmological constant are tuned to satisfy 
q . the Randall-Sundrum constraint on the brane. Braneworld models admit a 

c/2 . wider range of possibilities for dark energy than standard LCDM. In these 

models the luminosity distance can be both smaller and larger than the lu- 
minosity distance in LCDM. Whereas models with dh < (if, (LCDM) imply 
w = p/p > —1 and have frequently been discussed in the literature, mod- 
els with di > di (LCDM) have traditionally been ignored, perhaps because 
within the general-relativistic framework, the luminosity distance has this 
property only if the equation of state of matter is strongly negative (w < —1). 
Within the conventional framework, 'phantom energy' with w < — 1 is beset 
with a host of undesirable properties, which makes this model of dark en- 
ergy unattractive. Braneworld models, on the other hand, have the capacity 
to endow dark energy with exciting new possibilities (including w < —1) 
without suffering from the problems faced by phantom energy. For a sub- 



class of parameter values, braneworld dark energy and the acceleration of the 
universe are transient phenomena. In these models, the universe, after the 
current period of acceleration, re-enters the matter-dominated regime so that 
the deceleration parameter q(t) — ► 0.5 when t S> to, where to is the present 
epoch. Such models could help reconcile an accelerating universe with the 
requirements of string/ M-theory. 

I. INTRODUCTION 

One of the most sensational discoveries of the past decade is that our universe is cur- 
rently accelerating. This observation finds support in the luminosity measurements of high- 
redshift supernovae [1], measurements of degree-scale anisotropies in the cosmic microwave 
background [2] and, indirectly, in the observations of gravitational clustering [3] . Although a 
cosmological constant appears to satisfy all current observations, the formidable fine-tuning 
difficulties which accompany a non-evolving A-term have prompted theorists to investigate 
alternative, evolving forms of dark energy. Theoretical models of dark energy usually involve 
matter fields possessing unusual properties such as a negative equation of state. Popular 
dark-energy models include scalar-field-based 'quintessence' models, models based on quan- 
tum particle production, Chaplygin gas etc. (see [4-6] for recent reviews). In this paper, 
we examine a radically different mechanism for the late time acceleration of the universe in 
which dark energy effectively emerges from the gravity sector and not the matter sector of 
the theory. 

We study a higher- dimensional cosmology in which the observable universe is a four- 
dimensional 'brane' embedded in a five-dimensional bulk. Models of this kind appeared 
after the seminal paper by Randall & Sundrum (RS) [7]. Subsequent intensive investigation 
showed that cosmology based on the RS model exhibits departure from standard general- 
relativistic behaviour only at very early cosmological times [8-10]. The model which we study 
in this paper generalizes the RS model in that, besides the brane and bulk cosmological 
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constants, it also includes the scalar curvature term in the action for the brane. In an 
interesting recent development, Deffayet, Dvali, and Gabadadze (DDG) [11] demonstrated 
that the presence of the scalar curvature term in the action for the brane can lead to a late- 
time acceleration of the universe even in the absence of any material form of dark energy 
[11-13]. The main difference between our model and the DDG model is that, in addition 
to the scalar curvature term in the action for the brane, we also include the brane and bulk 
cosmological constants. (The bulk and brane cosmological constants were set equal to zero 
in the DDG model.) 

The presence of bulk and brane cosmological constants in our braneworld (in addition 
to the brane curvature term) leads to several qualitatively new features that distinguish our 
model from the DDG model as well as from the scalar-field-based 'quintessence' models. An 
important property of our braneworld universe is that the (effective) equation of state of 
dark energy can be w < — 1. (The DDG braneworld model and quintessence models have 
w > — 1.) In addition, the acceleration of the universe in our braneworld scenario can be 
a transient phenomenon. We find that a class of braneworld models accelerate during the 
present epoch but revert back to matter-dominated expansion at late times. A transiently 
accelerating braneworld does not possess an event horizon and could reconcile a currently 
accelerating universe with the demands of string/M-theory. 

Our paper is organised as follows. After describing the braneworld model in the next 
section, in Sec. Ill we proceed to the case of vacuum braneworlds, that is, solutions with zero 
stress-energy tensor of matter. We derive a general expression for the effective cosmological 
constant in such models and describe some classes of its symmetric solutions, including 
empty static homogeneous and isotropic universes. We also present simple expressions for 
a static universe filled with matter. In Sec. IV we study the cosmological evolution of the 
braneworld. We shall show that braneworld models have properties similar to the observed 
(accelerating) universe for a large region of parameter space. For instance, the model under 
consideration can mimic the evolution of a FRW universe in which dark energy has pressure- 
to-energy-density ratio pj p both larger as well as smaller than the critical value of —1. For a 
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subclass of parameter values, braneworld dark energy and the acceleration of the universe are 
transient phenomena so that the universe, after the current period of acceleration, re-enters 
the matter-dominated regime. In Sec. V we formulate our conclusions. 



II. BASIC EQUATIONS 

In this paper, we consider the case where a braneworld is the timelike boundary of a five- 
dimensional purely gravitational Lorentzian space (bulk), which is equivalent to the case of 
a brane embedded in the bulk with Z 2 symmetry of reflection with respect to the brane. 
The theory is described by the action [14,15] 

S = eM 3 \[ (K- 2A b ) -2 f k] + [ (m 2 R - 2a) + f L (h ab , 0) . (1) 

.J hulk J brane J J brane ^ ' J brane 

Here, 1Z is the scalar curvature of the metric g ab in the five-dimensional bulk, and R is 
the scalar curvature of the induced metric h ab = g a b — n a n b on the brane, where n a is the 
vector field of the inner unit normal to the brane, and the notation and conventions of 
[16] are used. The quantity K = K ab h ab is the trace of the symmetric tensor of extrinsic 
curvature K ab = /i c a V ' c n b of the brane. The symbol L(h ab , 0) denotes the Lagrangian density 
of the four-dimensional matter fields whose dynamics is restricted to the brane so that 
they interact only with the induced metric h ab . All integrations over the bulk and brane 
are taken with the natural volume elements \/ zr g d 5 x and \/ — h d 4 x, respectively, where g 
and h are the determinants of the matrices of components of the corresponding metrics in 
a coordinate basis. The symbols M and m denote, respectively, the five- dimensional and 
four-dimensional Planck masses, A b is the bulk cosmo logical constant, and a is the brane 
cosmo logical constant, also called the brane tension. 

The parameter e = ±1 in the first term of action (1) reflects the possibility of different 
relative signs between the bulk and brane terms in the action. In the context of homogeneous 
Friedmann cosmology that will be subject of this paper, the differential equations on the 
brane will not depend on the sign of e. However, different signs of e result in different ways 



4 



of embedding of a brane with one and the same induced metric in the bulk space, which 
may be important for studying perturbations of the theory. 

The term containing the scalar curvature of the induced metric on the brane with the 
coupling m 2 in action (1) is often neglected in the literature. However, this term is qualita- 
tively essential for describing the braneworld dynamics since it is inevitably generated as a 
quantum correction to the matter action in (1) — in the spirit of an idea that goes back to 
Sakharov [17] (see also [18]). Note that the effective action for the brane typically involves 
an infinite number of terms of higher order in curvature (this was pointed out in [14] for the 
case of braneworld theory; a similar situation in the context of the AdS/CFT correspon- 
dence is described in [19]). In this paper, we retain only the terms linear in curvature. The 
linear effects of the curvature term on the brane were studied in [12,21], where it was shown 
that it leads to four-dimensional law of gravity on sufficiently small scales. For some recent 
reviews of the results connected with the induced-gravity term in the brane action, one may 
look into [22]. 

Action (1) leads to the bulk being described by the usual Einstein equation with cosmo- 
logical constant: 

Gab + Kg ab = , (2) 

while the field equation on the brane is 

m 2 G ab + ah ab = r ab + eM 3 (K ab - h ab K) . (3) 

Here, r ab is the stress-energy tensor on the brane, and it stems from the last term in action 
(1). In Eq. (3) we have underlined the term whose presence makes braneworld theory 
different from general relativity. (It should perhaps be mentioned that the variation of 
action (1) which includes the Gibbons-Hawking surface term / K leads to the Israel 

J brane 

junction conditions [20] on the brane, as demonstrated in [14,15].) From (3) one can show 
that there exists a scale whose value determines the domain in which general relativity is 
approximately valid [12,21]. Indeed, we can write 
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G a6 ~rr 2 , K ab ~r 2 \ (4) 

where T\ and r 2 are the characteristic curvature radii of solution of Eq. (3). Then the ratio 
of the last term on the right-hand side of Eq. (3) to the first term on its left-hand side is of 
the order r\/ril, where i = 2m 2 /M 3 is a convenient definition to be used later. Thus, for 
r\ <C r 2 £ : the last term on the left-hand side of Eq. (3) can be neglected, and general relativity 
is recovered. In typical situations, e.g., in cosmology, we often have r\ ~ r 2 ~ r, and general 
relativity becomes valid for r <C £. This fact was used in [12], where the cosmological model 
with o" = and A b = was considered. Therefore, an important distinguishing feature of 
the braneworld models which we consider, is that they can show departure from standard 
FRW behaviour at late times and on large scales. This property should be contrasted with 
the braneworld cosmology based on the Randall-Sundrum model [7] (with m — 0), in which 
nonstandard behaviour is encountered at very early times [8-10]. 
By contracting the Gauss identity 

R abc d = hJh h vh c k h d {R Sgk i + K ac K b d - K bc K a d (5) 

on the brane and using Eq. (2), one obtains the 'constraint' equation 

R - 2A b + K ab K ab - K 2 = , (6) 

which, together with (3), implies the following closed scalar equation on the brane: 
M 6 (R - 2A b ) + (m 2 G ab + ah ab - r ab ) (m 2 G ah + ah ab - r ab ) - ^ (m 2 R - 4a + rf = , 

(7) 

where r = h ab r ab . One method for obtaining solutions of the theory consists in first solving 
the scalar equation (7) on the brane together with the stress-energy conservation equation, 
and then integrating the Einstein equations in the bulk with the given data on the brane 
[9,23]. 

We note that Eq. (7) describes the evolution on the brane in terms of its intrinsic quan- 
tities, and that all homogeneous and isotropic cosmological solutions for the brane can be 
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obtained from this equation [23], as will also be done below. All such solutions are embed- 
dable in the Schwarzschild-AdS five-dimensional bulk. However, one cannot confine oneself 
to this single equation in the general case, e.g., in studying perturbations of the cosmological 
solution, where one must also solve equations in the bulk imposing certain boundary and/or 
regularity conditions. From this general viewpoint, Eq. (7) cannot be regarded as complete: 
only those of its solutions are admissible which can be developed to regular solutions in 
the bulk. This illustrates the importance of the precise specification of boundary and/or 
regularity conditions in the bulk for a complete formulation of the brane- world theory [23]. 

The gravitational equations in the bulk can be integrated by using, for example, Gaussian 
normal coordinates, as described, e.g., in [24]. Specifically, in the Gaussian normal coordi- 
nates (x,y), where x = {x a } are the coordinates on the brane and y is the fifth coordinate 
in the bulk, the metric in the bulk is written as 

ds\ = dy 2 + h a/3 (x, y)dx a dx 13 . (8) 

Introducing also the tensor of extrinsic curvature K a b of every hypersurface y = const in the 
bulk, one can obtain the following system of differential equations for the components h a/3 
and K a p \ 

— - ^ = K* p - KK% - -6% (R + 2A b + K\K\ - K 2 ) 
dy 6 V J (g) 

= R a p-KK a p-p a p A h , 



dh a/3 



= 2h on ICp, (10) 



dy 

where R a p are the components of the Ricci tensor of the metric h a p induced on the hy- 
persurface y = const, R = R a a is its scalar curvature, and K = K a a is the trace of the 
tensor of extrinsic curvature. The second equality in (9) is true by virtue of the 'constraint' 
equation (6). Equations (9) and (10) together with the 'constraint' equation (6) represent 
the 4+1 splitting of the Einstein equations in the Gaussian normal coordinates. The initial 
conditions for these equations are defined on the brane through Eq. (3). We emphasise once 



again that, to obtain a complete braneworld theory in the general case (including a stability 
analysis), one must also specify additional conditions in the bulk such as the presence of 
other branes or certain regularity conditions. In our paper, we deal only with the homo- 
geneous and isotropic cosmology on the brane, so this issue does not arise. In this sense, 
we are studying here the cosmological features common to the whole class of braneworld 
models described by action (1) with arbitrary boundary conditions in the bulk. 

III. VACUUM BRANES AND STATIC BRANES 

We proceed to the cosmological implications of the braneworld theory under consider- 
ation. In this section, we discuss the situation pertaining to a vacuum brane, i.e., when 
the matter stress-energy tensor r ab = 0. It is interesting that the brane approaches this 
condition during the course of cosmological evolution provided it expands forever and its 
matter density asymptotically declines to zero. In this case, Eq. (7) takes the form 

{m % + Imj R + m 4 (R ab R ab - ^R 2 ^ - 4M 6 A RS = , (11) 

where 

Ars = - + ^- (12) 

Rb 2 3M 6 V ; 

It is important to note that the second term in Eq. (11) has precisely the form of one of the 
terms in the expression for the conformal anomaly, which describes the vacuum polarization 
at the one- loop level in curved space-time (see, e.g., [18]). 1 It therefore immediately follows 
that all symmetric spaces are solutions of Eq. (11) with appropriate A RS , just as they are 
solutions of the Einstein equations with one-loop quantum-gravitational corrections [25]. 



l It is interesting that, while the conformal anomaly term R a t,R ab — ^R 2 cannot be obtained by the 
variation of a local four-dimensional Lagrangian, the very same term is obtained via the variation 
of a local Lagrangian in the five-dimensional braneworld theory under investigation. 



8 



Symmetric spaces satisfy the condition V a Rbcde = 0, which implies that geometrical invari- 
ants such as R abcd R abcd , R ab R ab , and R are constants so that Eq. (11) becomes an algebraic 
equation. Prominent members of this family include: 

(i) The homogeneous and isotropic de Sitter space-time 



ds 2 = -dt 2 + — cosh 2 Ht 



d X 2 + sin 2 X (d0 2 + sin 2 Odcj) 2 ) 



(13) 



where — oo < t < oo, < x, < ir, < < 2tc. The four dimensional metric (13) has 
the property R a b = 3H 2 h a b , and formed the basis for Starobinsky's first inflationary model 
sustained by the quantum conformal anomaly [26]. 

(ii) The homogeneous and anisotropic Nariai metric [27,28] 



ds 2 = k 2 [-dt 2 + cosh 2 tdr 2 + d0 2 + sin 2 Odcf) , 



(14) 



where k = constant, — oo < t < oo, < 9 < n, < < 27r and for which R a b = h a b /k 2 . 
In fact, it is easy to show that any metric for which R and R a bR ab are constants will 
automatically be a solution to Eq. (11) with an appropriate choice of A RS - 

Both de Sitter space and the Nariai metric belong to the class of space-times which 
satisfy the vacuum Einstein equations with a cosmological constant 



R 



ab 



Ah 



ab ■ 



(15) 



Such space-times also satisfy Eq. (11) if 



A = — 

77T 



3M 6 



(2m 2 



/W6 \ 2 



(16) 



Equation (16) expresses the resulting cosmological constant on the brane in terms of the 
coupling constants of the theory. For the frequently discussed special case m — 0, one 
obtains A = Ars. The two signs in (16) correspond to the two different ways in which the 
lower-dimensional brane can form the boundary of the higher- dimensional bulk [14,13]. In 
the case of a spherically symmetrical bulk, the ' — ' sign corresponds to a brane as a boundary 
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for which the inner normal to the brane points in the direction of decreasing bulk radial 
coordinate. 

The condition Ars = is the well-known fine-tuning condition of Randall and Sundrum 
[7] and leads to the vanishing of the cosmological constant on an empty brane if we set 
m = in (1). Note that, under the Randall-Sundrum condition, expression (16) with the 
sign opposite to the sign of the quantity 3M 6 /2m 2 + a also gives a zero value for the resulting 
cosmological constant on the brane, but the other sign usually leads to A ^ 0. 

We would like to draw the reader's attention to the fact that Eq. (16) is meaningful only 
when the expression under the square root is nonnegative. When it is negative, solutions 
describing the corresponding empty universe simply do not exist. This leads to the following 
important conclusion: a universe which contains matter and satisfies 

3M 6 A RS , N 

— ?>1, 17 

(3M 6 /2m 2 + af V ; 

cannot expand forever. 

For the special case 3M 6 /2m 2 + a = 0, the expression for A is given by 

M 3 i 

A = ±— v/-3A RS . (18) 

In this case, both o and Ars must be negative in order that the corresponding empty universe 
exist, but the resulting cosmological constant on the brane can be of any sign. 

Another interesting example is that of a static empty universe. The radius (scale factor) 
a of such a universe is easily determined from (11) to be 

where k = ±1 is the sign of the spatial curvature. One can see that the radius of the 
universe can be arbitrarily large. In the general case, the development of this solution to 
the five-dimensional bulk leads to a Schwarzschild-anti-de Sitter metric. It was shown in 
[15] that, for k = 1, this metric is purely anti-de Sitter (with zero Schwarzschild mass) if 
the constants of the theory satisfy the condition 
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a A h 3M 



rn 



- + 7 = 0, (20) 

2 4m 4 v ; 



which implies negative brane tension a. It should be pointed out that the static and empty 
braneworld solution described by (19) does not possess a general-relativistic analog, since, 
in general relativity, a static cosmological model (the 'static Einstein universe') cannot be 
empty (see, for instance, [4]). Furthermore, from (19) we find that the static empty universe 
can be spatially open (k, = —1) — for example, in the case Ars < and o > — 3M 6 /2m 2 , 
— again a situation without an analog in general relativity. 

For static homogeneous and isotropic braneworlds filled with matter, Eq. (7) gives the 
following relation: 



a 2 



PtotUot + 3p to t) - 3A b M 6 = 3k m 2 ( Ptot + 3 Ptot ) - 3M 6 , (21) 



where the total energy density p tot and pressure p tot include the contribution from the brane 
tension, i.e., 

Ptot = p + a, p^t =p-a, (22) 

and k = 0, ±1 corresponds to the sign of the spatial curvature. This relation reduces to (19) 
for p = p = 0. 

Having obtained all these solutions on the brane, one can find the corresponding solutions 
in the bulk by integrating Eqs. (9) and (10) with the initial conditions on the brane given by 
Eq. (3), as described in Sec. II. In doing this, one can consider various additional conditions 
in the bulk, for example, existence of other branes, or one can impose certain regularity 
conditions. It is worth noting that one and the same cosmological solution on the given 
brane can correspond to different global solutions in the bulk, for example, other branes 
may be present or absent, be static or evolving, close or far away from our brane, etc. In the 
most general case (for instance in the absence of special symmetries on the brane) integration 
on the brane needs to be performed in conjunction with dynamical integration in the bulk. 
All such situations must be separately studied and issues such as their stability to linearized 
perturbations must be examined on a case-by-case basis. 
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IV. COSMOLOGICAL CONSEQUENCES OF BRANEWORLD DARK ENERGY 



A. Cosmological evolution 

In the homogeneous and isotropic cosmological setting, Eq. (7) can be integrated with 
respect to time with the result [23] 



K 



p + a 

3m 2 



(23) 



where C is the integration constant that corresponds to the black-hole mass of the 
Schwarzschild-(anti)-de Sitter solution in the bulk, H = a/a is the Hubble parameter on 
the brane, and k — 0, ±1 corresponds to the sign of the spatial curvature on the brane. A 
generalization of this equation to the case of absence of Z 2 symmetry was obtained in [29,23]. 
For the frequently considered special case m — 0, this equation reduces to the familiar one 
[8] 



H 2 + 



K 



Ab 
6 



C_ 



(p + af 



(24) 



a' b a* 9M 6 
Equation (23) can be solved with respect to the Hubble parameter with the result [14,13]: 

2 



H z + 



p + o 
3m 2 



+ 



1±\ l + £ 2 



p + cr At 



3m 2 



6 



C_ 

7a 



(25) 



or, equivalently, 



# 2 + 4 

a 1 



Ab C 1 
— - H 1 

6 a 4 £ 2 



p + cr 
3m 2 



(26) 



where the length scale 



1 = 



2m 2 
IP 



(27) 



was introduced earlier in Sec. II. Again, the 'i'signs in (25) and (26) correspond to two 
different ways of bounding the Schwarzschild-(anti)-de Sitter bulk space by the brane [14,13]: 
the lower sign corresponds to the case where the inner normal to the brane (which bounds 
the bulk) points in the direction of decreasing bulk radial coordinate. Alternatively, the 
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two different signs in (25) and (26) could correspond to the two possible signs of the five- 
dimensional Planck mass M. Henceforth we shall refer to models with the lower '— ' sign as 
BRANE1 and models with the upper '+' sign as BRANE2. 

In the ensuing discussion, we shall neglect the decaying term C/a 4 in Eq. (25) (called 
dark radiation), assuming it to be negligibly small at present. Similarly to [11], we also 
introduce the dimensionless cosmological parameters 

O p0 O — O - O _L O A_ 

m ~3mm^ K ~ a 2 H 2 ' a ~ 3m 2 H 2 ' k ~ PH* ' Ab ~ 6H 2 ' 1 8J 

where the subscript '0' refers to current values of cosmological quantities. The system of 
cosmological equations with the energy density p dominated by dust-like matter can now be 
written in the transparent form using Eq. (25): 

• BRANE1 

H 2 (z) 



H 2 



a n (i+^) 3 + n K (i+z) 2 + + 2n e - 2^n e ^/n m (i+z) 3 + n fT + n e + n Ah , 

(29) 



BRANE2 



= n m (i+z) 3 + n K (i+ z ) 2 + n a + 2n e + 2Jn e Jn m (i+ z y + n a + n e + n Ab , 

H 1 ? 

(30) 



where z is the cosmological redshift. We have underlined the terms which cause our equa- 
tions to differ from their general-relativistic counterparts. The constraint equations for the 
cosmological parameters are most conveniently derived using Eq. (26). We obtain 

i - n K + n Ab = (j/n e + n m + + n Ab ± y/n^j 2 . (31) 

This equation can be transformed to a more convenient form. In doing this, we note that 
the theory makes sense only in the case Vti + f2 m + il a + il Ab > 0. Since VLt > 0, when taking 
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the square root of both sides of this equation, we must consider two possibilities for the sign 
of the quantity f2 m + Q a + Q Ab under the square root in (31). 
In the case 

Q m + Q a + Q Ab > , (32) 
taking the square root of both sides of (31) and rearranging terms, we get 

- n K + n Ab =f = \M + fi m + n a + n Ab . (33) 

Now taking the square of both sides of the last equation, we finally obtain the constraint: 

fi m + n K + n a ± 2 V / q^ v / i - n K + n Ah = 1 . (34) 

In the opposite case 

-n t < Q m + Q a + Q Ah <0, (35) 
taking the square root of both sides of (31) and rearranging terms, we get 

y/i - n K + n Ab - Jn e = ±^/n e + n m + n a + n Ab . (36) 

Taking the square of both sides of the last equation, we finally obtain the constraint: 

n m + q k + n a + 2^n e ^/i -n K + n Ah = i. (37) 

To summarise, in the case Q m + Q a + Q Ab > 0, the constraint equations in the BRANE1 
and BRANE2 models read: 

ft m + n K + vt a - 2y/n t y/i -n K + n Ab = i (branei), (38) 

tt m + tt K + tt a + 2^n £ \/l-tt K + tt Ab = 1 . (BRANE2) (39) 

In the opposite case Q m + Q a + ^A b < 0, both constraint equations for BRANEI and 
BRANE2 are given by Eq. (39). Again, we have underlined the terms which cause our 
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equations to differ from their general-relativistic counterparts. It is easy to see that the 
general- relativistic limit of this theory is reached as il e — > (i.e., M — > 0). In this case, 
setting Q K = 0, Q a = Q\, we find that (29) & (30) give the expression for the Hubble 
parameter in the LCDM model: 



(40) 



An example of the behaviour of the Hubble parameter in the braneworld models is illustrated 
in Fig. 1. 



BRANE1 
BRANE2 




redshift z 



FIG. 1. The Hubble parameter in units of Hlcdm(z) is plotted as a function of redshift for the 
two braneworld models BRANE1 and BRANE2. Whereas H(z) in BRANE2 is larger than H(z) 
in LCDM, in the case of BRANE1 the value of H(z) is smaller than its value in LCDM. Parameter 
values are: Q K = 0, tt e = 1.0, Q m = 0.3, and Q Ah = 1, 10, 10 2 , 10 3 (top to bottom for BRANE2 
and bottom to top for BRANE1). The value of Q a is determined from (38) Sz (39). Parameter 
values for the LCDM model are Q m = 0.3, S7a = 0.7. For large values of the bulk cosmological 
constant ^A b , BRANE1 and BRANE2 become indistinguishable from LCDM. 

It is convenient to rewrite Eq. (12) in terms of the dimensionless cosmological parameters 

Arc fi? 



RS 



IRS 



4fi. 
15 



(41) 



from where it follows that the Randall-Sundrum constraint Ars = has the simple form 

^-4QAv b = 0. (42) 

This constraint implies f2A b > (A b < 0) and ensures the vanishing of the effective cos- 
mological constant in the future for the braneworld embedding described by the BRANE1 
model for Q a > or for Q a < and Vt t > Q Ah , and by the BRANE2 model for Q a < and 
Qe < Q\ h - The case where the Randall-Sundrum constraint is imposed will be discussed in 
detail in Sec. IV C below. It will be shown that, in the case of BRANE1 model, for the case 
Q K = 0, Eq. (38) would imply that fl m > 1, which apparently contradicts current observa- 
tions of the large-scale distribution of matter. This is an important result since it suggests 
that for reasonable values of the matter density (fi m <C 1) and Q K — 0, the BRANE1 model 
demands the presence of a nonvanishing effective cosmo logical constant. The value of the 
cosmological constant is given by Eq. (16) and can be determined by using observational 
data to fix the constants of the theory. On the other hand, the BRANE2 model is in prin- 
ciple observationally compatible with the Randall-Sundrum constraint, which can result in 
vanishing effective cosmological constant in the future thus making the current acceleration 
of the universe a transient phenomenon. 

Typical values of the f2 parameters (28) that we consider in this paper are of order 
fl ~ 1. For such values, the fundamental constants of the theory have the following orders 
of magnitude: 

m 2 ~ Ml ~ 10 19 GeV , M ~ 100 MeV , A b ~ ^- ~ Hi ~ 10~ 56 cm" 2 . (43) 

m 2 

The smallness of the bare cosmological constants in the bulk and on the brane represents a 
fine-tuning similar to what is the case for the cosmological constant in the standard LCDM 
model. However, even with such small values of the bare cosmological constants in action 
(1), the braneworld evolution exhibits some qualitatively new properties (discussed in the 
following subsections) when compared to the case where these constants are set to zero. 
This is the reason why we find it important to keep open the possibility of nonzero, albeit 
small, values of A b and a. 
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A possible procedure for testing the braneworld model against observations is presented 
in the appendix. 



B. Cosmological tests of braneworld models 



Observations of high-redshift type la supernovae indicate that these objects are fainter 
than they would be in a standard cold dark matter cosmology (SCDM) with H m = 1 [1]. 
This observation is taken as support for a universe which is accelerating, fuelled by a form of 
energy with negative pressure (dark energy). In standard FRW cosmology the acceleration 
of the universe is described by the equation 



where the summation is over all matter fields contributing to the dynamics of the universe. 
It is easy to show that a necessary (but not sufficient) condition for acceleration (a > 0) 
is that the strong energy condition is violated for at least one of the matter fields in (44), 
so that p + 3p < 0. In the case of the popular ACDM model this requirement is clearly 
satisfied since p m = in pressureless (cold) matter, while p^ = —p^ = —A/8nG in the 
cosmological constant. The situation with respect to braneworld models is different since, 
as we have demonstrated in the previous section, braneworld evolution is distinct from FRW 
evolution at late times. However it is easy to show that braneworld models can accelerate. 
We demonstrate this by noting that a completely general expression for the deceleration 
parameter q = —d/aH 2 is provided by 



where H(z) is given by (29)-(39) and the derivative is with respect to z. The current value 
of the deceleration parameter is easily calculated to be 



a 



(44) 



a 




(45) 




(46) 
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where the lower and upper signs correspond to BRANE1 and BRANE2 models, respectively. 
Thus the present universe will accelerate (q < 0) for brane parameter values which satisfy 



i"-( 1± V n. + iV?n« + nJ <1 - (4?) 

It should be pointed out that the proposition to use the induced curvature term on a 
braneworld to account for the accelerated phase of our universe was first made in [13] in the 
context of the braneworld model of [12], which is the special case of our BRANE2 model 
with il a = f2 Ab = 0. It was subsequently discussed and tested in [11,30]. In the present 
paper, we allow both Q a ^ and fljs, b 7^ 0, and the presence of these two free parameters 
makes the braneworld model more flexible and allows for qualitatively new behaviour which 
shall be discussed in detail in this section. The connection between large extra dimensions 
and dark energy has also been investigated in [31]. 

Observationally, a pivotal role in the case for an accelerating universe is played by the lu- 
minosity distance di{z), since the flux of light received from a distant source varies inversely 
to the square of the luminosity distance, T oc d~[ 2 . This effect is quantitatively described by 
the magnitude-luminosity relation: uib = M + 25 + 51og 10 <iL, where tub is the corrected 
apparent peak B magnitude and M is the absolute peak luminosity of the supernova. A 
supernova will therefore appear fainter in a universe which possesses a larger value of the 
luminosity distance to a given redshift. 

In a FRW universe, the luminosity distance is determined by the Hubble parameter and 
three-dimensional spatial curvature [4]: 

d L {z) = ^ + Z S (77o - 77) , (48) 

-^Oy I "total — 1| 

where 

dz' 



770 - 77 = tfoVl^totai - 1| j Q -^-yy , (49) 

and S(x) is defined as follows: S(x) — sinx if k — 1 (Jltotai > 1), S(x) = sinhx if k — — 1 
(^totai < 1), and S(x) = x if k = (f^otai = !)■ 
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Inflationary models suggest that the universe is nearly flat with Q K ~ and we shall 
work under this assumption in the rest of this paper (see also [32]). In this case, Eq. (48) 
simplifies to 

d L (z) dz' 

m = y ^M' (50) 

where H (z) is determined by (29) for BRANE1, and by (30) for BRANE2. In Fig. 2 
we show the luminosity distances for the BRANE1 & BRANE2 models. Also shown for 
comparison is the value of di(z) in a spatially-flat two-component FRW universe with the 
Hubble parameter 

H{z) = H [0,(1 + zf + n x (l + z) 3 ^] 1/2 , (51) 

where fl x describes dark energy with equation of state w = px/px- Three cosmo logical 
models will be of interest to us in connection with (51): 

(i) SCDM: The standard cold dark matter universe with Q m = 1 and Q x = 0. 

(ii) LCDM: Cold dark matter + a cosmological constant with w — — 1. 

(iii) Phantom models: Cold dark matter + 'phantom energy' satisfying w < — 1 [33]. 

We find from Fig. 2 that the luminosity distance in both braneworld models exceeds that 
in SCDM. In fact, BRANE1 models have the unusual feature that their luminosity distance 
can even exceed that in LCDM (for a fixed value of Q m ) ! In fact it can easily be shown that 

d L dS {z) > d L BRANE1 (z) > d L LCDM (z) , (52) 

where d^ s (z) refers to the luminosity distance in the spatially flat coordinatization of de Sit- 
ter space (equivalently, the steady state universe). The second inequality presumes a fixed 
value of Q m . BRANE2 models show complementary behaviour 

d L LCDM (z) > d L BKANE2 (z) > d L SCDM (z) , (53) 

where the first inequality is valid for a fixed value of Q m . In the case O = 0, the equa- 
tions of the braneworld theory formally reduce to those of general relativity, and we have 

^BRANEl^) = rfz BRANE2 (2) = ^LCDM^) 
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FIG. 2. The luminosity distance is shown as a function of redshift for the two braneworld models 
BRANE1 & BRANE2, LCDM, SCDM, and 'phantom energy'. All models, with the exception of 
SCDM, have Q m = 0.3. SCDM has ft m = 1. The BRANE1 & BRANE2 models have Q e = 0.3 and 
vanishing cosmological constant in the bulk. LCDM and the phantom model have the same dark 
energy density f^A = Qx = 0.7. The equation of state for dark energy is wa = —1 for LCDM and 
w = px/px = —1-5 for phantom. The luminosity distance is greatest for BRANE1 & phantom, 
and least for SCDM. BRANE1 & BRANE2 lie on either side of LCDM. 

One might add that the behaviour of BRANE1 is mimicked by FRW models with w < — 1, 
whereas BRANE2 resembles dark energy with — 1 < w < [4]. In fact, from Fig. 2 we see 
that the luminosity distance in the BRANE1 model is quite close to what one gets from 
'phantom energy' described by (51) with w = —1.5. (The parameters for this BRANE1 
model are Q m = Qe = 0.3, = 0, and Q a = 1 — Q m + 2\fTl^ « 1.8.) It should be pointed 
out that phantom energy models were introduced by Caldwell [33], who made the impor- 
tant observation that dark-energy with w < — 1 appears to give a better fit to supernova 
observations than LCDM (which has w = —1). Unfortunately, phantom models have sev- 
eral bizarre properties, which may be the reason why cosmologists have been reluctant to 
take these models seriously despite their seemingly better agreement with observations (see 
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however [34]). Some strange properties of phantom energy with a very negative equation of 
state (w < —1) are summarised below (see also [33,35]): 

(i) A negative equation of state suggests that the effective velocity of sound in the medium 



v = \J\dp/dp\ can become larger than the velocity of light. 

(ii) The expansion factor of a universe dominated by phantom energy grows as 

-, 2/3(1+™) 



a(t) ~ a (t eq ) 



[1 + w)— w 

teq 



, w<-l, (54) 



where t eq marks the epoch when the densities in matter and phantom energy are equal: 
Pm{t e q) — Px(t eq )- It immediately follows that the scale factor diverges in a finite amount 
of cosmic time 

a(t) ^00 as t — > t* — (^j^ — ) ^ et i ' (^) 

Substitution of z — > — 1 and w < — 1 in (51) shows that the Hubble parameter also diverges 
as t — > £*, implying that an infinitely rapid expansion rate for the universe has been reached 
in the finite future. 

As the universe expands, the density of phantom energy (w < —1) grows instead of 
decreasing (w > — 1) or remaining constant (w = —1), 



p(t) oc 



. t 

[l + w) w 

teq 



(56) 



reaching a singular value in a finite interval of time p(t) — > 00, t — > £*. This behaviour 
should be contrasted with the density of ordinary matter which drops to zero: p m — > 
universe dominated by phantom energy is thus doomed to expand towards a 
physical singularity which is reached in a finite amount of proper time. (An exact expression 
for the time of occupance of the phantom singularity is given in [36] , which also contains an 
interesting discussion of related issues.) 

At this stage one must emphasize that, although the BRANE1 model has several features 
in common with phantom energy, (which makes us believe that, like the latter, it too is 
likely to provide a good fit to supernova data), it is not necessarily afflicted with phantom's 
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FIG. 3. The distance modulus (m— M) of Type la supernovae (the primary fit of the Supernova 
cosmology project) is shown relative to an empty S7 m — > Milne universe (dashed line). The solid 
line refers to the distance modulus in BRANE1 with fi^ = fi m = 0.3, and vanishing cosmological 
constant in the bulk. The dotted line (below the solid) is LCDM with (f&A,O m ) = (0.7,0.3). The 
uppermost and lowermost (dot-dashed) lines correspond to de Sitter space (^A,O m ) = (1,0) and 
SCDM (fi A ,fim) = (0,1), respectively. 

pathologies. Indeed, in a broad range of parameters, both BRANE1 and BRANE2 are 
physically well motivated and remain well behaved during all times. The present paper 
therefore adds an important new dimension to the current debate about the acceleration of 
the universe by showing that cosmological models with cIl(z) > d^ CDM (z) are possible to 
construct within the framework of the braneworld scenario and should be taken seriously. 
Future studies will address the important quantitative issue of whether braneworld models 
of dark energy provide as good (or better) fit to high-z supernova observations as LCDM (see 
Fig. 3 as an illustration). (It is appropriate to mention here that the weak energy condition 
p + p > can also be violated in a class of scalar-tensor theories of gravity, as discussed in 



We should add that the 'angular-size distance' is related to the luminosity distance 



[37].) 
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FIG. 4. The age of the universe (in units of the inverse Hubble parameter) is plotted as a 
function of the cosmological redshift for the models discussed in Fig. 2. (The phantom model is 
not shown.) BRANE1 models have the oldest age while SCDM is youngest. 

d L through d A = d L (l + z)~ 2 , therefore much of the above analysis carries over when one 
discusses properties of the angular-size distance within the framework of braneworld models. 
Some cosmological features of braneworld models are shown in Figs. 4-6. In Fig. 4, the age 
of the universe at a given cosmological redshift 



is shown for the two braneworld models and for LCDM & SCDM. We find that the age of 
the universe in BRANE1 (BRANE2) is larger (smaller) than in LCDM for identical values 
of the cosmological density parameter f2 m . This is a direct consequence of the fact that 
the Hubble parameter in BRANE1 (BRANE2) is smaller (larger) than in LCDM. Both 
braneworld models are significantly older than SCDM. 

Considerable insight into the dynamics of the universe is provided by the cosmological 
deceleration parameter (45). Our results, shown in Fig. 5, indicate that at late times the 
BRANE1 (BRANE2) universe accelerates at a faster (slower) rate than LCDM (with iden- 
tical fl m ). Curiously, the BRANE1 universe shows an earlier transition from deceleration to 




(57) 
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FIG. 5. The deceleration parameter q(z) is shown for BRANE1, BRANE2 and LCDM. The 
model parameters are as in Fig. 2. For reference it should be noted that q = 0.5 for SCDM while 
de Sitter space has q = — 1. 

acceleration than any of the other models. (For the given choice of parameters this transition 
takes place at z ~ 1 for BRANE1 and at z ~ 0.7 for LCDM. The BRANE2 model begins 
accelerating near the present epoch at z ~ 0.) A related point of interest is that at z > 2 
the deceleration parameter in BRANE1 marginally exceeds that in SCDM indicating that 
the BRANE1 model is decelerating at a faster rate than SCDM (q = 0.5). In conventional 
models of dark matter this behaviour can occur only if the equation of state of the dark 
component is stiffer than dust, implying w > in (51), or if the universe is spatially closed. 
On the other hand, the current acceleration rate of BRANE1 in our example (q ~ —0.7) 
significantly exceeds that of LCDM (q ~ —0.55) with an identical value of Q m = 0.3 in 
both models. Within the framework of four-dimensional Einstein gravity, this situation can 
only arise if the equation of state of dark energy is strongly negative: w < — 1 in (51). 

The unusual high-z behaviour of the deceleration parameter in BRANE1 can be better 
understood if we consider the cosmological density parameter 
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FIG. 6. The dimensionless matter density £l m (z) is shown for the two braneworld models and 
LCDM. (f2 m = 1 in SCDM.) Parameter values are the same as in previous figures. BRANE1 has 
the interesting feature that £l m (z) exceeds unity for z > 1. 



where H(z) is given by (29)-(39) for braneworld models. From Fig. 6 we notice that, for 
z>\, the value of Q m (z) in BRANE1 exceeds its value in SCDM (f2 m = 1). This is precisely 
the redshift range during which q , (^)branei > q(z)scr>M- Thus, the rapid deceleration of 
BRANE1 at high redshifts can be partly attributed to the larger value of the matter density 
Q m (z) at those redshifts, relative to SCDM. 

Having established partial similarity of BRANE1 with phantom models at low redshifts, 
we can investigate the analogy further and calculate the effective equation of state of dark 
energy 



where £l m (z) is given by (58). One notes that w(z) has a pole-like singularity at z ~ 1 
for BRANE1, which arises because Q m (z) crosses the value of unity at z ~ 1 (see Fig. 6). 



Ho 



(58) 



H(z) 




(59) 
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This demonstrates that the notion of 'effective equation of state' is of limited utility for 
this model. Equations (29), (30), (45), and (59) also illustrate the important fact that dark 
energy in braneworld models, though similar to phantom energy in some respects, differ 
from it in others. For instance, in both braneworld models, w(z) — > —0.5 at z ^> 1 and 
w(z) — > — 1 as z — > — 1, whereas phantom energy has w(z) < —1 at all times. 

A useful quantity is the current value of the effective equation of state of dark energy in 
braneworld theories: 



where the lower and upper signs, as usual, correspond to BRANE1 and BRANE2 models, 
respectively. We easily see that w < — 1 for BRANE1, whereas w > — 1 for BRANE2. 



The braneworld models under consideration also admit the intriguing possibility that 
the current acceleration of the universe may not be a lasting feature. It may be recalled 
that most models of dark energy including the cosmological constant have the property that, 
once the universe begins to accelerate, it accelerates forever. As shown in a number of recent 
papers, an eternally accelerating universe is endowed with a cosmological event horizon 
which prevents the construction of a conventional S-matrix describing particle interactions 
within the framework of string or M-theory [39]. In this section we show that, provided 
the Randall-Sundrum constraint relation (42) is satisfied, the acceleration of the universe 
can be a transient phenomenon in braneworld models. An anisotropic solution of Bianchi V 
class with the same feature was described in [40]. 

From Eqs. (29) and (30) we obtain the following asymptotic expressions for the Hubble 
parameter Hoo as z — > — 1, assuming that the universe expands forever: 




(60) 



C. Disappearing Dark Energy 




(61) 
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FIG. 7. The deceleration parameter is plotted as a function of redshift for the BRANE2 model 
with the Randall-Sundrum constraint (42), fl Ab = 2, and (Q m , Sl t ) = (0.3, 1.2), (0.2, 1.6), (0.1, 1.98) 
(top to bottom). The vertical (dot-dashed) line at z = marks the present epoch, while the 
horizontal (dashed) line at q = corresponds to a Milne universe [a(t) oc t] which neither accelerates 
nor decelerates. Note that the universe ceases to accelerate and becomes matter dominated in the 
future. 

where the lower and upper signs correspond to BRANE1 and BRANE2 models, respectively. 
In applying the Randall-Sundrum constraint (42), we first consider the case where Q a > 0. 
Then 



and 



-"0 



(62) 



(63) 



One can see that this expression vanishes for the lower sign. Thus, for positive Q a , it is 
the BRANE1 model that leads to vanishing effective cosmological constant in the future. 
However, in this case, the constraint equation (38) with Q K = becomes 



2 v /^( v /i + ^ Ab - yjn^) = 1 
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FIG. 8. The effective equation of state for dark energy in the BRANE2 model is shown as 
a function of redshift. Model parameters are as in the previous figure. Note that the past and 
future asymptotes of w(z) are quite different: w(z) — ► — 1/2 for z » 0, while w(z) — > for 
z — > — 1. Braneworld dark energy therefore effectively disappears in the future, giving rise to a 
matter-dominated universe. 

and implies f2 m > 1, which is hardly compatible with the observations. 
In the case of Vl a < 0, we have 



(65) 



and 



IT x 2 

-"0 



(66) 



If > f2A b j then this expression vanishes for the lower sign, which brings us back to the 
nonphysical BRANE1 models with f2 m > I. The case tie > 1 + fU b is compatible with 
the condition Q m + Q a + ^A b < 0, which corresponds to BRANE1 (29) with the constraint 
equation (39). However, it can be shown that these conditions imply Q m < 0, which is also 
physically unacceptable. 
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There remains the case of Vt a < and Vt t < {l\ h . In this case, expression (66) vanishes 
for the upper sign, which corresponds to BRANE2 models. The constraint equation (39) 
with Q K = now reads 



and implies f2 m < 1. 

Therefore, BRANE2 with Q a < and il e < f2 Ab , provides us with an interesting ex- 
ample of a physically meaningful cosmological model in which the current acceleration of 
the universe is a transient phenomenon. An example of this behaviour as probed by the 
deceleration parameter is shown in Fig. 7, which demonstrates that the current period of 
cosmic acceleration takes place between two matter-dominated epochs. We emphasize that 
these models require negative brane tension a. Since an observer in this model resides on 
a negative-tension brane one must ponder over the issue of whether such a braneworld will 
be perturbatively stable and hence physically viable. We consider this to be an open ques- 
tion for future investigations. Remarks made at the end of Sec. Ill are relevant, however, 
since one and the same cosmological solution on the 'visible' (negative tension) brane can 
correspond to many different global conditions in the bulk, for instance, other ('hidden') 
branes may be present or absent, static or evolving, close to or far away from our brane, 
etc. Moreover, here the sign of the parameter e = ±1 in action (1) is expected to be of 
crucial importance, since, in particular, the bulk gravity will be localised around the neg- 
ative tension (visible) brane in the case of e = — 1. Issues of perturbative stability involve 
the examination of all such distinct global solutions on a case-by-case basis. Such a study, 
though interesting and important, clearly lies beyond the scope of the present paper. 

Useful insight into the BRANE2 model is also provided by the effective equation of state 
of dark energy (59). Our results shown in Fig. 8 indicate that the past and future behaviour 
of dark energy in the braneworld universe can be very different. The past behaviour w(z) — > 
—0.5 for z ^> 1 arises because, in a spatially flat braneworld, the second most important 
contribution to braneworld expansion at high redshifts is caused by the (1 + z) 3 ^ 2 term in 




(67) 
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(30); see also [11]. The future behaviour w(z) — > as z — > —1, on the other hand, reflects 
the decreasing importance of dark energy as the universe expands. The acceleration of the 
universe is therefore a transient phenomenon which ends once the universe settles back into 
the matter-dominated regime. 

Finally, we should mention that a transiently accelerating regime also arises in a class of 
BRANE2 models which do not satisfy the Randall-Sundrum constraint (42). In these models 
the current epoch of acceleration is succeeded by an epoch during which the deceleration 
parameter grows without bound. This unusual 'future singularity' is reached in a finite 
interval of expansion time and is characterised by the fact that both the matter density and 
the Hubble parameter remain finite, while a — > oo (a feature that distinguishes it from the 
phantom singularities discussed in Sec. IV B). A detailed discussion of the 'new' singularities 
which occur in braneworld models can be found in [41]. 

V. CONCLUSIONS 

In this paper, we have considered cosmological implications of a braneworld model de- 
scribed by action (1), which contains both bulk and brane curvature terms and cosmological 
constants. The curvature term for the brane arises naturally, as a quantum correction 
from the matter part of the brane action, and significantly changes the behaviour of the 
braneworld theory. For example, as is well known, braneworld cosmology without this term 
deviates from general relativity at large matter densities, i.e., at early cosmological times; 
see Eq. (24). In the model that we considered, on the other hand, early cosmological evo- 
lution remains virtually unaffected for a broad range of parameters, significant deviations 
from standard cosmology appearing only during later times. 

We restrict our attention to the important case where the brane forms the boundary 
of the five-dimensional bulk space. This is equivalent to endowing the bulk with the Z 2 
reflection isometry with respect to the brane. The presence of the curvature term in the 
action leads to two families of braneworld models. These two families (called BRANE1 & 
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BRANE2) differ in the manner in which the brane forms the boundary of the five- dimensional 
bulk. For example, if a spatially three-spherical brane is embedded into the three-spherically 
symmetric bulk, then one can discard either the interior part of the bulk, or its exterior 
(asymptotically flat) part. In this example, the models which discard the exterior part 
leaving the interior were classified as BRANE1 models, while the complementary models 
were called BRANE2. Similar classification arises in the case of spatially flat and open 
braneworld models. Alternatively, the two different families of braneworld models can be 
regarded as corresponding to the two possible signs of the five- dimensional Planck mass M. 

Braneworld models of dark energy considered in this paper have the interesting and un- 
usual property that their luminosity distance di can exceed that in LCDM . This is unusual 
since, within the general-relativistic framework, the luminosity distance has this property 
only if the equation of state of dark energy is strongly negative (w < —1). Caldwell [33] 
has shown that dark energy with w < — 1 (phantom energy) may in fact fit supernova 
observations better than LCDM. However, phantom energy is beset with a host of un- 
desirable properties which makes this model of dark energy unattractive. We show that 
braneworld models have all the advantages and none of the disadvantages of phantom mod- 
els and therefore endow dark energy with exciting new possibilities. A recent analysis of 
braneworld models in [38] has demonstrated that BRANE1 models (which generically have 
w < — 1) provide an excellent fit to supernovae observations for higher values of the matter 
density (Q m > 0.3), whereas lower values (Q m < 0.25) are preferred by BRANE2 models 
(which generically have w > —1). 

A distinctive feature of the braneworld scenario discussed in this paper is that it al- 
lows for a universe which is transiently accelerating. Recent investigations indicate that an 
eternally accelerating universe, which possesses a cosmological event horizon, prevents the 
construction of a conventional S-matrix describing particle interactions within the frame- 
work of string or M-theory [39]. We have demonstrated that braneworld models can enter 
into a regime of accelerated expansion at late times even if the brane tension and the bulk 
cosmological constant are tuned to satisfy the Randall-Sundrum constraint on the brane. In 
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this case, braneworld dark energy and the acceleration of the universe are transient phenom- 
ena. In this class of models, the universe, after the current period of acceleration, re-enters 
the matter-dominated regime. We have shown that viable models realising this behaviour 
are those of BRANE2 type. (Since these braneworlds have a negative brane tension the 
question of their stability is important and needs to be examined in detail. We shall return 
to this issue in a companion paper.) Thus braneworld models can give rise to a transiently 
accelerating phase thereby reconciling a dark energy dominated universe with the require- 
ments of string/M-theory. (A similar observation was previously made in [40] in the context 
of an anisotropic cosmological solution of Bianchi V class.) 
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APPENDIX: TESTING THE MODEL 

In this appendix, we present a possible procedure for testing the braneworld model 
against observations. 

There are the following four possibilities to be considered: 

(i) n m + n a + n Ab >o [E q . (32)]. 

(a) yjl - Q K + Q Ah - y/n~ £ > =>• {(29), (38)} and {(30), (39)} are possible. 



(b) yjl -n K + Vl Ah - VTT e < =>• only {(29), (38)} is possible, 
(ii) Q m + Q a + Q Ab < [Eq. (35)]. 

(a) yJl-£l K + tt Ah - ^/U~ e > =>• only {(30), (39)} is possible. 



(b) yjl - Q K + Q Ab - < =>• only {(29), (39)} is possible. 
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To see these possibilities, one must consider Eqs. (33) and (36). 

Let us introduce the quantity F = ^Jl — Q K + Q Ab —y/Qe that enters the above conditions 
(a) and (b). Now we note that, when using (38) for determining fi CT , we have 

n m + n a + n Ab = i-n K + n Ab + 2^n e ^Ji-n K + n Ab , (ai) 

so that condition (i) is satisfied automatically for any sign of F. Similarly, when using (39) 
for determining Q a , we have 

n m + n a + n Ab = i-n K + n Ab - 2y/n ey /i -n K + n Ab , (A2) 

so that condition (ii) is satisfied automatically for F < 0. This means that conditions (i) 
and (ii) need not be checked in the procedure. 

The case of the Randall-Sundrum constraint with Q K = is much simpler, since it allows 
only BRANE2 model. Note that the condition il e < £l Ab must be satisfied in this case. [See 
the explanations above Eq. (67).] The value of £lg can be obtained from (67): 

2y/^= 1-an ■ (A3) 

y 1 + ^A b - y"A b 

So that 

Sl t <n Ab =s> 1 - ft m + 2Q Ah < + Q Ab . (A4) 

This inequality can be simplified by taking the square of both sides: 

(i - n m + 2n Ab ) 2 < 4^ Ab (i + n Ah ) (i - n m ) 2 - 4Q m Q Ab < o . (as) 

Finally, we get the inequality 

< A6 > 

Summarizing these results, one can suggest the following procedure for testing the 
braneworld models: 
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Specify the values of fie , fl\ b , fl K , and fl m (initial fTs). In doing this, note that always 
fie > by definition, < fl m < 1 from observations, and 1 — fl K + fi Ab > by physical 
consistency [see Eq. (31)]. In addition, one can assume the universe to be spatially flat: 
fl K = 0. 

(A) To study the BRANE1 model (29) with constraint (38): 
Use (38) to calculate fl a , and (29) to calculate di(z). 

(B) To study the BRANE1 model (29) with constraint (39): 

In specifying the values for the initial fTs, always choose Q e > 1 — Q K + fl Ah . Use (39) 
to calculate fl a , and (29) to calculate di(z). 

(C) To study the BRANE2 model (30) with constraint (39): 

In specifying the values for the initial f2's, always choose fie < 1 — f2 K + f2 Ab . Use (39) 
to calculate Q a , and (30) to calculate cIl(z). 

(D) The case of disappearing dark energy discussed in Sec. IV C and described by BRANE2 
model is studied separately as follows: 

Set f2 K = 0, specify the value of f2 m in the range < Q m < 1, and choose f2A b in the 
range given by (A6). Now use (67) to find fie , and then (65) to find fl a . 

Finally, use (30) to calculate d L (z). 
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